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New results in heat-kernel asymptotics on manifolds with boundary 
1. Introduction 



The investigation of the heat equation for differential (or pseudo-differential) operators on 
manifolds with boundary remains of high relevance in physics and mathematics. Mathe- 
maticians are more interested in the general properties of spectral geometry and spectral 
asymptotics, with the aim of interpreting the various heat-kernel coefficients with the help 
of invariance theory [1-3] in a neat, elegant and deep way, and of finding resolvent and 
heat kernel also when the boundary- value problem is pseudo-differential [4] . Heat-equation 
proofs of the index theorem are also available, by now [3,5]. Theoretical physicists, on the 
other hand, remain more interested in the applications, e.g. boundary conditions for 
Euclidean quantum gravity [6-13], one-loop semiclassical approximation [8,12], and the 
quantization programme for gauge theories on manifolds with boundary [14-17]. 

In the first part of our presentation, we shall consider an m-dimensional Riemannian 
manifold, say (M,g), a vector bundle V over M, with a connection V, and operators of 
Laplace type: 

P = -g ab V a V b - E, (1.1) 

with E an endomorphism of V. The corresponding heat kernel is, by definition, a solution, 
for t > 0, of the equation 

(j t +P^U(x,x';t) = 0, (1.2) 
which obeys the initial condition 

lim / U (x, x'\ t)p(x') y/det g dx = p(x), (1.3) 
' '"./.\/ 

and suitable (local or non-local) boundary conditions 

BU(x,x';t)} =0 . (1.4) 

J dM 

In the case of local boundary conditions, we shall assume them to be of the form [17] 

n o \( [ip] dM 



a i-n; \[v NV ] dM 1 "°' (L5) 
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where IT is a self-adjoint projection operator, A is a tangential differential operator on the 
boundary of M : 



A= (i-n) 



^(rv, • v,r) + s 



(i-n), (1.6) 



and cp are the fields, i.e. the smooth sections of the bundle V. With the notation of Eq. 
(1.6), V is the induced connection on <9M, P are endomorphism- valued vector fields on 
the boundary, and S is an endomorphism of the vector bundle over dM which is a copy 
of |V] 9M ; with sections given by [<p]dM- and S are anti-self-adjoint and self-adjoint, 
respectively, and are annihilated by IT on the left and on the right, i.e. TIT 1 = F l U = 
IIS = SU = 0. As is shown in [3,11,12,16,17], one arrives at such boundary conditions 
whenever one tries to obtain gauge- and BRST-invariant boundary conditions in quantum 
field theory. 

We will be interested in the asymptotic expansion as t — > + of the L 2 -trace [20] 

Tr L 2(fe- tp ) = / f(x)Tr v U(x,x;t)y/det]jdx, (1.7) 

where / G C°°(M). This is also called the integrated heat kernel. Equation (1.7) means 
that one first takes the fibre trace of the heat-kernel diagonal. Composition with the 
smearing function /, and integration over M, yields the functional trace (1.7). The results 
for the original boundary- value problem are eventually recovered by setting / = 1, but it 
is crucial to keep the smearing function arbitrary throughout the whole set of calculations. 
In other words, we consider the global, integrated asymptotics, for which 

oo 

Tr L2 (/e"* p ) ~(47rt)"* jN*i4 5 (/,P,B). (1.8) 

The coefficients A™ can be always expressed in the form 

A 9 (f, P, B) = C s (/, P) + B 9 (f, P, B). (1.9) 

The coefficients C» and Bq are the interior part and the boundary part, respectively. The 
interior part vanishes for all odd values of n, whereas the boundary part only vanishes if 
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n = 0. The interior part is obtained by integrating over M a linear combination of local 
invariants of the appropriate dimension, built universally and polynomially from the met- 
ric, the Riemann curvature R a bcd of M, the bundle curvature, say Q a &, the endomorphism 
E in the operator (1.1), and their covariant derivatives. By virtue of the Weyl theorem on 
the invariants of the orthogonal group [1,2], which is 0{m) for M, these polynomials can 
be found by using only tensor products and contraction of tensor arguments. Moreover, 
the boundary part is obtained upon integration over dM of another linear combination of 
local invariants. In that case, however, the structure group is 0(m— 1), and the coefficients 
of linear combination are universal functions, independent of m, unaffected by conformal 
rescalings of the metric g, and invariant in form (i.e. they are functions of position on the 
boundary, whose form is independent of the boundary being curved or totally geodesic). 

Thus, the general form of the A a coefficient is a well posed problem in invariance the- 
ory, where one has to take all possible local invariants built from /, R a bcd , O a &, Kij, E, S, F l 
and their covariant derivatives (hereafter, is the extrinsic-curvature tensor of the 
boundary), integrating eventually their linear combinations over M and <9M, respectively. 
For example, in the boundary part B™(f,P,B), the local invariants integrated over dM 
are of dimension n — 1 in tensors of the same dimension of the second fundamental form 
of the boundary, for all n > 1. The universal functions associated to all such invariants 
can be found by using functorial methods, e.g. conformal variations, lemmas on product 
manifolds [1,2], jointly with the consideration of some particular manifolds (see below). 

Section 2 shows how to apply the conformal-variation method to heat-kernel asymp- 
totics for the generalized boundary-value problem resulting from (1.1), (1.5) and (1.6). 
Section 3 describes the recent results on the issue of strong ellipticity for this class of local 
boundary- value problems. Fourth-order operators on manifolds with boundary are studied 
in section 4, and some recent ideas on integro-differential boundary conditions in Euclidean 
quantum gravity are discussed in section 5. Section 6 is devoted to some open problems, 
and a number of new heat-kernel coefficients are evaluated in the appendix. 
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2. Conformal variations and heat-kernel coefficients 



The conformal- variation method [18-21] is based on the behaviour of the heat kernel under 
conformal rescalings. The idea is to perform a conformal deformation of the differential 
operator P and of the boundary operator £>, with a deformation parameter, say e (see 
below). This is used to determine a set of recurrence relations (see (2.1)-(2.3)) which 
reduce the evaluation of heat-kernel asymptotics to the solution of a system of algebraic 
equations for a finite set of coefficients. More precisely, the "conformal variations" of the 
metric and of the inward-pointing normal are g a b(s) = e 2e ^ g a b(0), N a (e) = e~ £ f N a (0). 
The conformal variation of the potential terms E(e),S(e) and P(e) is chosen in such a 
way that the operator P and the boundary operator B transform according to P(e) = 
e- 2sf P(0),B(e) = e~ sf B(0). This implies, in particular, that P(e) = e" e/ r i (0). Three 
basic conformal- variation formulae hold [1-3]: 



d_ 

de 



A 9 (l, e- 2 ^P(0), e-"B(0))l = (m - n)A 9 (f, P(0), B(0)) , (2.1) 

J s=0 



^4 S (l, P(0) - eH, e" e 'B(0)) J = ^9-1 P ^), B(0)), (2.2) 



de 2 



if,e- 2e/ P(0),e- £/ i3(0)) 



= 0, 



£=0 



(2.3) 



where H is another function e C°°(M), different from /. Our analysis, relying on Refs. 
[20,21], is restricted to the case when the following conditions hold: 



o, 



(2.4) 



(2.5) 



ViP = 0. 



(2.6) 



In other words, the matrices P commute with each other, the matrix T 2 = TjP commutes 
with S, and the matrices T l are taken to be covariantly constant with respect to the induced 
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connection, V, on the boundary. These assumptions are indeed quite restrictive, but, for 
the time being, not even the A\ coefficient is known if (2.4)-(2.6) fail to hold, because 
formidable technical difficulties are found to arise [17,20,21]. Thus, Eqs. (2.4)-(2.6) do 
not express a compulsory step, but represent the current limits of what one is able to do. 

Although the Eqs. (2.4)-(2.6) impose severe restrictions, a non-trivial invariance 
theory is found to contribute many new local invariants to the heat-kernel asymptotics. 
Explicitly, one finds the following structure of the first few interior terms [18,20,21]: 



Bi(/,P,B) = / Tr( 7 /), 

JdM 



(2.7) 



B 1 (f,P,B) = ^ f Tr\f(b K + b 2 S) + b 1 f. !N + f(a 1 K ij T i rA 
O JdM L v 7 



(2- 



Bs (/, P, B) = ^ J Tr {/ (c E + Cl R + c 2 R m lN + c 3 K 

+ f ;N (c 5 K + CgS) + C 6 / ;ArA r| 



+c 4 KijK ij +c 7 SK + c 8 S 2 



384 J dM 



+ \iK i jT i T j s + N jN rv j + n 2 R l aj rT J ' 



+hn iN r 



+ f3 1 f. N K ij T i T^ 



(2.9) 



For any subsequent interior term, the number of new local invariants resulting from the 
occurrence of T l in the boundary operator is higher and higher (but finite). For example, in 
the integrand for the coefficient B 2 (f, P, B), 33 new local invariants (involving F l ) multiply 
/, 7 local invariants multiply f ; N, and 1 local invariant, KijY l Y\ multiplies f-NN [20,21]. 
All universal functions occurring in (2.7)-(2.9) have been evaluated, thanks to the work in 
Refs. [15,17,18,20,21]. They are all expressed in terms of the functions 



Vl + T 2 , V /Z f2, Arctanh( V /Z T2). 
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For example, one finds, with the help of Eq. (2.1), and considering a particular case (i.e. 
a flat background with a totally geodesic boundary) which provides further information 



A remarkable piece of work, in Ref. [21], has evaluated all universal functions in (2.9), by 
using the formulae (2.1)-(2.3), jointly with a lemma on product manifolds [1,2], and the 
consideration of particular manifolds, i.e. the bounded generalized cone, and the manifold 
with B 2 x T m ~ 2 topology (T m ~ 2 being the equilateral (m — 2)-dimensional torus). The 
consideration of particular cases is always useful and, indeed, strictly necessary (so far), 
by virtue of the universal nature of the functions of T 2 one is looking for (see comments 
after (1.9)). 

3. Strong ellipticity 

If one studies compact Riemannian manifolds without boundary, it is enough to make 
sure that the leading symbol of the differential operator under consideration is elliptic. 
In the presence of boundaries, however, one has also to check that the strong ellipticity 
condition holds. At a technical level, this requires that a unique solution should exist of 
the eigenvalue equation for the leading symbol cr^P), subject to a decay condition at 
infinity and to suitable boundary conditions. To obtain a precise formulation, one has 
to consider local coordinates on the boundary <9M, say x k (k = 1, m — 1), the normal 
geodesic distance to <9M, say r, cotangent vectors on the boundary, say Q £ T*(dM) 
(j = 1, ...,m — 1), a real parameter, say u>, the graded leading symbol of the boundary 
operator in Eq. (1.5) (cf. the first matrix therein): 



[18,20]: 




(2.10) 




(3.1) 



and an arbitrary pair of boundary data, i.e. 




(3.2) 
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where such a ip' is, strictly, a smooth section of an "auxiliary" vector bundle over dM, 
endowed with a decomposition into sub-bundles, with half the dimension of the bundle of 
boundary data. 

By definition, strong ellipticity holds with respect to the cone C — 7Z + if a unique 
solution exists, say (p, of the equation [2,17,22] 

a L (P- x\ r = 0, -id r ) ~ Al] <p(r) = 0, (3.3) 
subject to the asymptotic condition 

lim ip(r) = 0, (3.4) 
and to the boundary condition (hereafter, x = {x } X = {Cj}) 

VC e T*(dM), VA G C-TZ+, V(C, A) ^ (0, 0) and V?//. What is crucial, for the physicists who 
are interested in the applications to quantum field theory (as well as for mathematicians 
who might be more interested in heat- kernel theory), is that the lack of strong ellipticity 
leads to a fibre trace of the heat-kernel diagonal which acquires a non-integrable part near 
the boundary. It is hence impossible to make sense of the integrated heat kernel (cf. (1.7)) 

Tr L2 (e"* p )= / Tr y t/(x,x;t) v / det g dx, (3.6) 

with the corresponding global asymptotics as t — > + (unless one studies a smeared form 
along the lines of (1.7)). 

This is why, in Ref. [17], a systematic investigation of strong ellipticity for local 
boundary- value problems involving operators of Laplace and Dirac type has been carried 
out. For the former class of differential operators, which are our main source of concern 
so far, the work in Refs. [17,22] may be summarized as follows. 
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(i) Let P be an operator of Laplace type, with boundary operator as in Eq. (1.5): 



Bp = 



n 
A i-n 



The boundary- value problem (P, Bp) is strongly elliptic with respect to C — 7Z + if and only 



if, for all £ 0, the matrix — iT^Q is positive-definite, where \(\ = + A/CjC J - 

(ii) The boundary-value problem for Euclidean Yang-Mills theory at one-loop level, with 
the gauge-invariant boundary conditions (N being the normal to the boundary) 



(5 a b - N a N b ) <p b 



dM 



= o, 



(3.7) 



(3.8) 



VV a = 0, 

L J dM 

is strongly elliptic with respect to C — 1Z + . 

(iii) In Euclidean quantum gravity at one-loop level, the vector bundle V is the bundle of 
symmetric rank- two tensor fields ip a b over M, with fibre metric 



E ab cd _ g a(c g d)b + ag ab g cd^ 



(3.9) 



where a is a real parameter different from — — (for a = — — , no inverse of E can be 
defined). On considering the tensor q a b = g a b — N a Nb, the projection operator which is 
self-adjoint with respect to the bundle metric (3.9) reads 



TT cd —— 

LL ab — 



Q(a Qb) 



a 



(a + 1) 



N a N b q 



cd 



(3.10) 



The boundary conditions invariant under infinitesimal diffeomorphisms on metric pertur- 
bations can be therefore expressed in the form [17] 



ncd , „ 
ab 'fed 



dM 



0, 



E abcd (a) V b <p cd 
9 



= 0. 



dM 



(3.11) 
(3.12) 
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If a = — |, the operator on metric perturbations is then of Laplace type, say P again. Re- 
grettably, one then finds that the boundary-value problem (P, -Bp), with Bp the boundary 
operator of the type (1.5) and (1.6) which gives rise to (3.11) and (3.12) with a = — |, 
is not strongly elliptic with respect to C — 7Z + . The work in Ref. [17] has also evaluated 
the non-integrable contribution to the fibre trace of the heat-kernel diagonal, which can 
be expected on general ground if strong ellipticity is violated, as we said after Eq. (3.5). 

In other words, only for Euclidean quantum gravity based on the Einstein action the 
strong ellipticity condition is violated, on using gauge-invariant boundary conditions of the 
form (3.11) and (3.12). On the other hand, as shown in Ref. [16], such boundary conditions 
can be re-derived, independently, by requiring invariance under BRST transformations. It 
therefore seems that a BRST-invariant quantization of the gravitational field presents 
serious inconsistencies on manifolds with boundary, unless one accepts the view according 
to which the BRST invariance of the amplitudes does not force the boundary conditions 
to be, themselves, BRST invariant [14]. 

4. Fourth-order operators on manifolds with boundary 

The current attempts to develop quantum field theories of fundamental interactions have 
led to the consideration of fourth-order or even higher-order differential operators on closed 
Riemannian manifolds [23-27], or on manifolds with boundary [3,28,29]. The analysis of 
the transformation properties under conformal rescalings of the background metric, say g, 
leads, in particular, to the consideration of conformally covariant operators, say Q, which 
transform according to the law 

Q u = e -( m + 4 ^/ 2 Q(uJ = 0)e( m - 4)w / 2 , (4.1) 

if g rescales as g w = e 2uj g, m being the dimension of the Riemannian manifold which is 
studied. One of the physical motivations for this analysis lies in the possibility to use the 
Green functions of such operators to build the effective action in curved space-times [25]. 

Another enlightening example is provided by the ghost sector of Euclidean Maxwell 
theory in vacuum in four dimensions. The corresponding field equations are well known 
to be invariant under conformal rescalings of g. On the other hand, the supplementary 
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(or gauge) conditions usually considered in the literature are not invariant under confor- 
mal rescalings of g. Even just in flat Euclidean four-space, conformal invariance of the 
supplementary condition is only achieved on making the Eastwood-Singer choice [30]: 

V b V b V c A c = 0, (4.2) 

where A c is the electromagnetic potential (a connection one-form in geometric language). 
The preservation of Eq. (4.2) under gauge transformations of A c : 

f Ac = Ac + Vc /, (4.3) 

is achieved provided that / obeys the fourth-order equation 

□ 2 / = 0, (4.4) 

where Q 2 is the box operator composed with itself: 

□ 2 = v a v a v & v 6 . 

In the corresponding quantum theory via path integrals, one thus deals with two inde- 
pendent ghost fields (frequently referred to as the ghost and the anti-ghost), both ruled 
by O 2 , which is a fourth-order elliptic operator, and subject to the following boundary 
conditions: 

[e] d M = 0, (4.5) 

V N e] =0. (4.6) 

Remarkably, since one now deals with a fourth-order elliptic operator, it is insufficient to 
impose just Dirichlet or Neumann (or Robin) boundary conditions. One needs instead 
both (4.5) and (4.6), which are obtained from the following requirements: 

(i) Gauge invariance of the boundary conditions on [3,28]. 

(ii) Conformal invariance of the whole set of boundary conditions. 

(iii) Self-adjointness of the Q 2 operator (see below). 

11 



New results in heat-kernel asymptotics on manifolds with boundary 

Although it remains extremely difficult to build a consistent quantization scheme via path- 
integral formalism for the full Maxwell field in the Eastwood-Singer gauge (the gauge-field 
operator on Af, perturbations being, then, of sixth order [3,28]), the investigation of the 
ghost sector remains of considerable interest in this case. There is in fact, on the one 
hand, the need to understand how to quantize a gauge theory in a way which preserves 
conformal invariance at all stages (as we just said), and on the other hand the attempt 
to extend the recent work on conformally covariant operators [23-27] to the more realistic 
case of manifolds with boundary. 

For simplicity, we consider the squared Laplace operator acting on scalar fields on a 
flat Euclidean background, in the case when curvature effects result from the boundary 
only. Moreover, motivated by quantum cosmology and Euclidean quantum gravity, the 
boundary is assumed to be a three-sphere of radius a, say, or a pair of concentric three- 
spheres [8,12]. The former case, in particular, may be viewed as the limiting case when 
the wave function of the universe is studied at small three-geometries (i.e. as a — > 0), as 
shown in Ref. [31]. 

In our problem it is hence possible to expand the scalar field on a family of three- 
spheres centred on the origin, according to the familiar relation [32] 

CO 

e(x,r) = J2^n(r)Q {n) (x), (4.7) 

n=l 

where r G [0, a], are the scalar harmonics on a unit three-sphere, S 3 , and x are local 
coordinates on S s . Thus, one is eventually led to study a one-dimensional differential op- 
erator of fourth order, and this makes it clear why all the essential information is obtained 
by the analysis of the operator B = on a closed interval of the real line, say [0, 1]. The 
operator B is required to act on functions which are at least of class C 4 (see (4.18)), and 
the following definition of scalar product (anti-linear in the first argument) is considered: 

(u,v)= / u*(x)v(x)dx. (4.8) 
Jo 
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We now want to study under which conditions the operator B is self-adjoint, which means 
that it should be symmetric, and its domain should coincide with the domain of the adjoint, 
say BK For this purpose, we first study the relation between the scalar products (Bu,v) 
and (u,Bv). We have then to integrate repeatedly by parts, using the Leibniz rule to 
express 



d ( d 3 u* \ d ( d 2 u* dv\ d ( du* d 



dx \ dx 3 J ' dx \ dx 2 dx J ' dx \ dx dx 2 J ' dx \ dx 3 J 



d ( „ d 3 \ 



u 



This leads to 



(Bu,v) = 



d 3 



u 



n 1 



dx 3 



d 2 u* dv 
dx 2 dx 



n 1 



+ 



du* d 2 v 
dx dx 2 



n 1 



u 



dh_ 

dx 3 



n 1 



+ (u,Bv). 



(4.9) 



Bearing in mind that the adjoint, B\ of is again the operator it is thus clear that 
the condition (Bu,v) = (u,B'v) is fulfilled provided that both u G D(B) and v G D(B^) 
obey one of the following four sets of boundary conditions: 



(i) First option: 



u(0) = u(l) = u'(0) = u'(l) = 0, 
u(0) = u(l) = v'(0) = v'(l) = 0. 



(4.10) 
(4.11) 



(ii) Second option: 



u(0) = u(l) = u"(0) = u"{l) = 0, 
u(0) = v(l) = v"(0) = v"(l) = 0. 



(4.12) 
(4.13) 



(iii) Third option: 



u '(0) = u'(l) = u"'(0) = u'"{\) = 0, 
v '(0) = v '(i) = ^"(0) = = 0. 



(4.14) 
(4.15) 



(iv) Fourth option: 



u"(0) = = u w (0) = u w (l) = 0, 



(4.16) 
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v"(0) = v"(l) = v"'(0) = v"'(l) = 0. (4.17) 

In other words, if the conditions (4.10) and (4.11), or (4.12) and (4.13), or (4.14) and 
(4.15), or (4.16) and (4.17) are satisfied, the domains of B and of its adjoint turn out to 
coincide [29]: 

D(B) = D(B ] ) = {u:ue AC 4 [0, 1], (4.10) or (4.12) 

or (4.14) or (4.16) holds} . (4.18) 

With our notation, AC 4 [0, 1] is the set of functions in L 2 [0, 1] whose weak derivatives up 
to third order are absolutely continuous in [0,1], which ensures that the weak derivatives, 
up to fourth order, are Lebesgue summable in [0,1], and that all u in the domain are of 
class C 4 on [0,1]. Of course, symmetry of B is also obtained with the boundary conditions 
just described. 

In other words, four sets of boundary conditions, (i) or (ii) or (iii) or (iv), can be chosen 

,4 

to ensure self-adjointness of the operator Hereafter, we first consider the option (i), 
since it agrees with boundary conditions motivated by the request of gauge invariance and 
conformal invariance, if the scalar field is viewed as one of the two ghost fields of Euclidean 
Maxwell theory in the Eastwood-Singer gauge. We also stress again that nothing is lost on 

,4 

studying just the "prototype" operator . The one-dimensional fourth-order operator 
may take a more complicated form in some set of local coordinates (see (4.19)), but is 
always reducible to the operator on the real line (more precisely, a closed interval of 
3? in our problems). 

The definition and evaluation of functional determinants remains a topic of crucial 
importance in quantum field theory. Here the task is even more interesting, because we 
are studying a fourth-order elliptic operator on a manifold with boundary. As shown in 
Refs. [3,28], the resulting eigenvalue equation for the modes occurring in the expansion 
(4.7) turns out to be, on the Euclidean four-ball, 



d 4 6 d s (2n 2 - 5) d 2 (2n 2 + 1) d (n 2 - l) 2 



dr 4 t dr 3 r 2 dr 2 r 3 dr r 4 



e n = A n e n . (4.19) 
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Thus, on setting M = the solution of Eq. (4.19) is expressed by a linear combination 
of Bessel functions and modified Bessel functions [29], i.e. 

rM-A In{MT) + A Kn{MT) + A UMt) 4- 4 M^l) f4 201 

£n(Tj — n h A 2 , n h ^3,n h ^4,n • l 4 -^ U J 

r r t t 

Since the Euclidean four-ball consists of a portion of flat Euclidean four-space bounded by 
a three-sphere, the coefficients A2, n and A^ n have to vanish Vn > 1, to ensure regularity 
of e n at the origin. One is thus left with scalar modes of the form 

In (Mr) Jn(Mr) 
e n (r)=Ai )n hA 3i „ . (4.21) 

T T 

These massless modes are subject to the boundary conditions (see (4.5), (4.6) and (4.10), 
(4.11)) 

[Sn]dM = 0, (4.22) 
[de n /dr] dM = 0. (4.23) 

The Eqs. (4.21)-(4.23) lead to the eigenvalue condition (denoting by a the radius of the 
three-sphere) 

detf In{Ma) UMa) Vo (4 24) 

\-I n (Ma) + MaI' n {Ma) -J n (Ma) + MaJ' n (Ma) ) ' K > 

which guarantees that non-trivial solutions exist for the coefficients A\^ n and A^^ n in (4.21). 
As proved in Ref. [29], Eq. (4.24) leads to the following ((0) value: 

«0) = ~ (425) 

This is of some interest, because the one- loop analysis remains crucial in the course of 
studying quantum theory as a theory of small disturbances [33] of the underlying classical 
theory. 

In the two-boundary problem one studies instead a portion of flat Euclidean four-space 
bounded by two concentric three-spheres. This case is very interesting because it is more 
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directly related to the familiar framework in quantum field theory, where one normally 
assigns boundary data on two three-surfaces (it should be stressed, however, that unlike 
scattering problems we are considering a path-integral representation of amplitudes in 
a finite region). On denoting by a and 6, with a > 6, the radii of the two concentric 
three-sphere boundaries, we can consider the complete form (4.20) of our scalar modes, 
because no singularity at the origin occurs in the two-boundary problem, and hence all 
linearly independent integrals are regular, for all r G [6, a] . We now impose the boundary 
conditions (4.5) and (4.6), which lead to the eigenvalue condition 



det 



/ I n {Mb) K n (Mb) J n (Mb) N n (Mb) \ 

F In (Mb) F Kn (Mb) F Jn (Mb) F Nn (Mb) 

I n (Ma) K n (Ma) J n (Ma) N n (Ma) 

\F In (Ma) F Kn (Ma) F Jn (Ma) F Nn (Ma) ) 



(4.26) 



where, for Z = i", K, J or N, we define [29] 

F Zn (Mx) = -Z n (Mx) + MxZ' n (Mx). (4.27) 

The work in Ref. [29] proves that Eq. (4.26) leads to a vanishing C(0) value: 

C(0) = 0. (4.28) 

The result (4.28) is found to hold for all boundary conditions described in Eqs. (4.10)- 
(4.17). 

To sum up, the original contribution of Ref. [29] is as follows. 

(i) The boundary conditions for which the squared Laplace operator is self-adjoint have 

,4 

been derived for the first time, taking as prototype the operator on a closed interval of 
the real line. Interestingly, four sets of boundary conditions are then found to arise, and the 
option described by (4.10) and (4.11) coincides, if the field in (4.7) were a ghost field, with 
the boundary conditions obtained from the request of gauge invariance of the boundary 
conditions on A^, when the Eastwood-Singer supplementary condition is imposed. 

(ii) Given the fourth-order eigenvalue equation (4.19), the contribution of the correspond- 
ing eigenmodes to the one-loop divergence has been derived for the first time on the 
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Euclidean four-ball (see (4.25)), or on the portion of flat Euclidean four-space bounded by 
two concentric three-spheres (see (4.28)). 

In our opinion, the result (i) is crucial because no complete prescription for the quantization 
is obtained unless suitable sets of boundary conditions are imposed, and Eqs. (4.10)-(4.17) 
represent a non-trivial step in this direction. The result (ii) is instead relevant for the anal- 
ysis of one-loop semiclassical effects in quantum field theory. In other words, if one has to 
come to terms with higher order differential operators in the quantization of gauge theories 
and gravitation, it appears necessary to develop techniques for a systematic investigation 
of one-loop ultraviolet divergences, as a first step towards a thorough understanding of 
their perturbative properties. 

Some outstanding problems are now in sight. First, it appears interesting to ex- 
tend our mode-by-mode analysis to curved backgrounds with boundary. In this case, 
the fourth-order conformally covariant differential operator is more complicated than the 
squared Laplace operator, and involves also the Ricci curvature and the scalar curvature of 
the background. Second, one should use Weyl's theorem on the invariants of the orthogo- 
nal group to understand the general structure of heat-kernel asymptotics for fourth-order 
differential operators on manifolds with boundary. A naturally occurring question within 
that framework is, to what extent functorial methods can then be used to compute all 
heat-kernel coefficients for a given form of the differential operator and of the boundary 
operator (cf. Refs. [20,21]). Third, the recently considered effect of tangential derivatives 
in the boundary operator [15-18,20-22] might give rise to generalized boundary conditions 
for conformally covariant operators. The appropriate mathematical theory is still lacking 
in the literature, but would be of much help for the current attempts to understand the 
formulation of quantum field theories on manifolds with boundary. 

5. Non-local boundary conditions in Euclidean quantum gravity 

The last decade of efforts on the problem of boundary conditions in (one-loop) Euclidean 
quantum gravity has focused on a local formulation, by trying to satisfy the following 
requirements: 
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(i) Local nature of the boundary operators [6-9,11,12,16,17,20,22]. 

(ii) Operator on metric perturbations, say P, and ghost operator, say Q, of Laplace type 
[17]. 

(hi) Symmetry, and, possibly, (essential) self-adjointness of the differential operators P 
and Q [11,12]. 

(iv) Strong ellipticity of the boundary-value problems obtained from the operators P and 
Q, with local boundary operators B\ and B 2: respectively [17,22]. 

(v) Gauge- and BRST-invariance of the boundary conditions and/or of the out-in (one- 
loop) amplitude [6,7,9,11,12,14,16,17,22]. 

At about the same time, in the applications to quantum field theory and quantum gravity, 
non-local boundary conditions had been studied mainly for operators of Dirac type (see, 
however, Ref. [10]), relying on the early work by Atiyah, Patodi and Singer on spectral 
asymmetry and Riemannian geometry [34]. What is non-local, within that framework, is 
the separation of the spectrum of a first-order elliptic operator (the Dirac operator on the 
boundary) into its positive and negative parts. This leads, in turn, to an unambiguous 
identification of positive- and negative-frequency modes of the (massive or massless) Dirac 
field, and half of them are set to zero on the bounding surface [3,8,35]. 

On the other hand, non-local boundary conditions for operators of Laplace type had 
already been studied quite intensively in the literature, from at least two points of view: 

(i) The rich mathematical theory of pseudo-differential boundary-value problems, where 
both the differential operator P and the boundary operator B may be replaced by integro- 
differential operators [4]. 

(ii) Bose-Einstein condensation models, where integro-differential boundary operators lead 
to the existence of bulk and surface states [36] . 
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For example, if P is an operator of Laplace type, mathematicians have derived many 
properties of the boundary-value problem [4] 



Pu = f in 0, 

Tu = ip at T = <90, 
where the boundary operator T can take the form 

Tu = 7o-u + 7>, 

or, instead, 



Tu = 71U + S ~f u + T[u. 



With this notation, one has [4] 



1j u = 



{-id n ) J u 



..7 = 0,1, 



(5.1) 
(5.2) 

(5.3) 
(5.4) 

(5.5) 



where d n is the inward-pointing normal derivative. Moreover, Tq and T[ are integral 
operators going from O to <9fi, and the map Sq acts on functions on dfl. 

In the case of the gravitational field, inspired by Eqs. (5.1)-(5.5), we consider a scheme 
where the differential operator on metric perturbations remains of Laplace type (as well 
as the ghost operator), whereas the boundary conditions are of integro-differential nature. 
This means that the full boundary operator, say B ab cd , may be expressed as the sum of 
a local operator, say B ab cd , obtained from projectors and first-order differential operators 
(see (1.5) and (1.6)), and an integral operator going from the background four-manifold, 
M, to its boundary dM, so that the boundary conditions read 



B ab cd h cd (x) 



J dM 



B ab cd h cd (x) 



J dM 



+ 



T ab cd {x,x')h cd {x')d^' 



dM 



(5.6a) 



dM 



This notation is a bit too general. We may decide, following DeWitt [37] , that unprimed 
lower-case indices refer to the point x and primed lower-case indices refer to the point x' . 
This leads to 



~B ab cd h cd (x) 




~B ab cd h cd (x) 


+ 


7 




dM 




dM 


JdM 



Tab d hc'd'dH' 



(5.66) 



dM 
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which is the form of the boundary conditions chosen hereafter [13]. 

Since we are concerned, for simplicity, with operators of Laplace type in a flat four- 
dimensional background (all curvature effects result then from the boundary only), it is 
very important for us to understand the effect of integro-differential boundary conditions 
on such a class of operators. For this purpose, following Ref. [4], we remark that, after 
integration by parts, one finds the Green formula (unlike Ref. [4], we define the Laplace 
operator with a negative sign in front of all second derivatives), for P = A, u E D{P), and 
v in the domain D(P*) of the adjoint of P: 

{Pu, v)q A u, v^j =(u, Avj ^+(l(pu, pvj^, (5.7) 

where the Green matrix reads, in our case [4] 

W = z(j Q, (5.8) 

whilst p is the (Cauchy) boundary operator, whose action reduces to 

pu = (70 u,7iu). (5.9) 

The same property (5.9) holds for v G D(P*). Suppose now that the boundary conditions 
are expressed in the integro-differential form (5.3): 

7o« + TqU = at T. 
The term (lApu, pv^j in Eq. (5.7), which is equal to 

(upu,pvj = i(7iii,7o«)r + i(7ow,7i^)r, (5.10a) 
can be then re-expressed as 

(l4pu,pv^ = z(7i«,7ov)r + i(-TqU,^iv)t, (5.106) 
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which implies that P* , the (formal) adjoint of P, can be obtained by adding to A a singular 
Green operator, i.e. 

P*v = Av + zTo*7iv, (5.11) 
supplemented by the local boundary condition 

7ot> = at T. (5.12) 

By contrast, if the boundary conditions (5.4) are imposed: 

7iu + S 7o« + T[u = at T, 

which modify the standard Neumann case, it is convenient to re-express 71 u, at the bound- 
ary, in the form 

71 it = -S ^ u - T[u, (5.13) 
and insert Eq. (5.13) into Eq. (5.10a). This implies that the adjoint of P now reads 

P*v = Av + iT[*7 v, (5.14) 

subject to the local boundary condition 

7it> = at T. (5.15) 

In other words, we are discovering a property which is known to some mathematicians, 
but not so familiar to physicists: if an elliptic differential operator (here taken to be 
of Laplace type) is studied with integro-differential boundary conditions, its adjoint is a 
pseudo-differential operator, subject to local boundary conditions. 

Self-adjointness problems are properly formulated by studying the so-called realization 
of the operator P [4]. In our case, this means adding to the Laplacian a singular Green 
operator, and considering a trace operator which expresses the integro-differential boundary 
conditions. More precisely, a Dirichlet-type realization of P = A is the operator 

B d =(a+G d ) , (5.16) 

V /To 
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where 

Gd^K^ + G', (5.17) 

To = 7 o + T^. (5.18) 

The Ki operators, for i = 0,1, are completely determined by the requirement of self- 
adjointness. In technical language, they are called Poisson operators [4]. Indeed, the 
domains of B D and its adjoint coincide if and only if [4] 

K 1 = iTg, (5.19) 

G' = G'*. (5.20) 
Moreover, a Neumann-type realization of P = A is the operator 

B n =(a+G n ) , (5.21) 

where 

G N = K 0l0 + F', (5.22) 

Ti =7i + ^o7o + T{. (5.23) 
The domains of and its adjoint are then found to coincide if and only if [4] 

K Q = iT[\ (5.24) 

So = -S *, (5.25) 

F' = F'\ (5.26) 

In the case of the gravitational field, our boundary operator (5.66) corresponds to the 
integro-differential trace operator (5.18). The local boundary operator B ab cd is taken to 
be the one for which the following conditions are imposed on metric perturbations on a 
three-sphere boundary of radius a [7]: 



hij 



= 0, (5.27) 

dM 
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h 



Oi 



dM 



o, 



dh 00 6 d j- ' 



0, 



(5.28) 
(5.29) 



dM 



where r G [0, a]. Equations (5.27)-(5.29) express, to our knowledge, the only set of local 
boundary conditions which are of Dirichlet type on hij and hoi, and for which strong 
ellipticity of the boundary- value problem is not violated (cf. section 3). Since we only 
want to modify the Dirichlet sector of such boundary conditions [13], which is expressed 
by (5.27) and (5.28), we have to require that (see (5.66)) 



T 00 cd =0 VcV. 
Thus, we eventually consider the operator 

A = 

where S is of the type (5.18) in its ij and Oi components, i.e. 



(5.30) 



A + G 
S 



(5.31) 



S^hcdix) 



dM 



hij(x) 



dM 



+ 



f T ii dd! h dd ,dY! 

JdM 



dM 



Soi Cd h c d(x) 




hoi 


+ 


7 




dM 




dM 


JdM 



Toi C d hc'd'dT,' 



(5.32) 



(5.33) 



dM 



and of the type (5.29) (cf. Eq. (5.23)) in its normal component /loo, i-e. 
S 00 cd h cd (x) 



dM 



dh 00 6 d y 

-7T- + - h oo ~ -7r(g J hij) 
or t or 



(5.34) 



dM 



Moreover, A is the standard Laplacian on metric perturbations in flat Euclidean four- 
space, and G may be viewed as the direct sum of Go and Gn (cf. example 1.6.16 in Ref. 
[4]), with 

G D = Km + G", (5.35) 



Gn — F ', 
23 



(5.36) 
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subject to the self-adjointness conditions 



{Ki) l3 cd h cd 



J dM 



dM 



Tj,j c d h c 'd'dT,' 



(Ki) 0l cd h cd 



— i 



dM 



(Ki) 



dM 



Toi C d K'd'dT,' 



dM 



dM 



oo d h c d 



dM 



o, 



cd 



ab 



n , cd\ 
U ab J > 

F' ab cd =(F' ab cd )\ 



(5.37) 

(5.38) 

(5.39) 
(5.40) 
(5.41) 



Note that, by virtue of Eq. (5.30), the counterpart of T{ vanishes in our problem, and 
hence no counterpart of Eq. (5.24) has to be imposed. Further to this, the condition (5.25) 
is satisfied by virtue of the boundary condition (5.34). 

The contribution of Ref. [13], summarized in the present section, consists of the 
proposal that the boundary conditions (5.66) should be considered as a serious candi- 
date for non-local boundary conditions in Euclidean quantum gravity; moreover, the self- 
adjointness conditions (5.37)-(5.41) have been derived, inspired by a careful analysis of 
the results first derived in Ref. [4]. At least four outstanding problems are now in sight: 

(i) Can one build explicitly a class of bulk and surface states in Euclidean quantum gravity 
with non-local boundary conditions, inspired by the work in Ref. [36]? The idea is then to 
obtain mode-by-mode solutions of the eigenvalue equations for metric perturbations, and 
insert them into Eq. (5.66) for a given form of T ab c d . One then looks for solutions which 
decay rapidly away from the boundary (the surface states), or remain non- negligible (the 
bulk states). 

(ii) Can one study heat-kernel asymptotics with non-local boundary conditions for the 
gravitational field? 
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(iii) Can one perform a path-integral quantization, if the non-local choice (5.66) is made 
for the boundary data? What form of boundary conditions should be imposed on ghost 
fields? 

(iv) Can one obtain a non-local formulation of Euclidean quantum gravity, where both 
the differential operators and the boundary operators are replaced by suitable classes of 
pseudo-differential operators, taking full advantage of the recent progress in the field of 
functional calculus for pseudo-differential boundary- value problems [4]? 

6. Concluding remarks 

Our presentation, although incomplete for obvious length limits, has tried to open the 
window on what we regard as several promising research lines for the years to come. They 
can be summarized by a number of questions, as follows. 

(i) When tangential derivatives occur in the boundary operator, can one compute further 
heat-kernel coefficients for operators of Laplace type in the Abelian (see (2.4) and (2.5)) 
and non-Abelian cases? 

(ii) What about the general structure of heat-kernel asymptotics for non-minimal operators 
[2,38,39] and conformally covariant operators [23,24] on manifolds with boundary? 

(iii) Can one obtain a non-local formulation of Euclidean quantum gravity with the help of 
functional calculus for pseudo-differential boundary- value problems, along the lines of Refs. 
[4] and [13]? Is there a corresponding heat- kernel asymptotics? Can one build explicitly 
bulk and surface states in Euclidean quantum gravity? 

(iv) What lesson should one learn from the recent proof that there is lack of strong ellip- 
ticity in Euclidean quantum gravity, upon choosing completely gauge-invariant boundary 
conditions with gauge- field and ghost operators of Laplace type [17,22]? Can one then 
define a smeared functional trace (cf. (1.7)) which is not ad hocl Is it correct to conclude 
that BRST-invariant boundary conditions for the quantized gravitational field turn out to 
be incompatible with the need for a well defined elliptic theory? 



25 



New results in heat-kernel asymptotics on manifolds with boundary 

The first two problems that we have selected might seem very technical, but their 
solution would lead to a substancial advancement of knowledge in spectral geometry and 
in the general theory of operators on manifolds. The remaining problems lie, instead, 
at the very heart of any attempt to get a deeper understanding of Euclidean quantum 
gravity. The whole programme might seem as either too ambitious or too abstract to the 
pragmatic reader, but such a stand would be most unfortunate. In our opinion, there are, 
by contrast, several good reasons for continuing advanced research along these lines, i.e. 

(i) Quantum gravity via Euclidean path integrals makes it possible to deal with a framework 
where, quite naturally, partition functions are defined. This is, in turn, crucial if one wants 
to combine quantum theory, whose predictions are of statistical nature [3], with general 
relativity. 

(ii) Non-local properties are frequently met in the course of studying modern quantum 
field theories via path-integral or canonical methods (cf. Ref. [40]). Thus, a non-local 
approach to Euclidean quantum gravity appears essential to complete the current efforts. 

(iii) Spectral geometry plays a crucial role if one tries to get a thorough understanding 
of the one-loop semiclassical approximation [3,8,12]. In a space-time approach [41-44], 
such an approximation provides the "bridge" in between the classical world and the (as 
yet unknown) full quantum theory (via path integrals). At some stage, any alternative 
approach to quantum gravity should be able to make contact with what one knows from 
a perturbative evaluation of transition amplitudes [12]. 

(iv) So much has been learned from the heat-kernel approach to index theory and to the 
theory of eigenvalues in Riemannian geometry [2-4,12] that any new result may have a non- 
trivial impact on Euclidean quantum gravity. Some readers might feel unhappy with the 
idea of having to formulate the quantum theory in a Euclidean framework, where no time- 
evolution exists, the differential operators are elliptic and the geometries are Riemannian 
(rather than Lorentzian as in general relativity). However, if one is interested in the 
most fundamental structures, it is not bad news to realize that a framework exists where 
some problems become well posed. From this point of view, one has to go ahead as 
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long as possible with Euclidean problems, before being sure that quantum theory needs 
a Lorentzian framework. In other words, if a manifold picture retains an important role 
in quantum field theory, the problem is open of whether one should regard the elliptic 
boundary-value problems as the most fundamental tool. In the years to come, hopefully, 
further progress in the topics discussed so far might lead to a deeper vision. 

Appendix 

The consideration of special cases remains of considerable help in the investigation of heat- 
kernel asymptotics, because the local invariants in the integrand for heat-kernel coefficients 
are multiplied by universal functions, whose form remains the same for all smooth bound- 
aries (see section 1). Thus, if one is able to evaluate a subset of the general set of universal 
functions, one can use this result, jointly with the algebraic equations resulting from (2.1)- 
(2.3) and from the application of lemmas on product manifolds [1,2], until one eventually 
gets enough algebraic equations to evaluate all universal functions corresponding to a given 
heat-kernel coefficient [21]. 

Thus, following Ref. [20], we consider the case when all curvatures vanish: R a bcd = 
0, Q ao = 0, Kij = 0, whilst E, S and T l are covariantly constant, in that 



VE = 0, VS = 0, ViF = 0. 



(A.l) 



The coefficient A 5 / 2 is a purely boundary term, and from Ref. [20] one finds in our case 




(A.2) 



where is the normal derivative of / of order r (e.g., = J-nn), and p n ^ 
functions generated by the following algorithm: 



are universal 




u 



k 



[u 2 + u + 



W + r 2 )]^ 



(A.3) 




(A.4) 
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1 



Pn,k 



2na n -i j k- n — T a n+ i } k- 



(A.5) 



n\(k — n)\ 

All universal functions are well defined and analytic provided that 1 + r 2 > 0. As remarked 
in Refs. [18,20] this reflects a crucial property: when T 2 = —1, the strong ellipticity of the 
boundary-value problem no longer holds. In our paper, we always assume that 1 + T 2 is 
positive. 

Now the explicit form of A 5 / 2 is, from Eq. (A. 2), 



Idi 



l 5/2 = / T^v 
>dM 



f( P 4,4S 4 + p 2 , 2 ES 2 + ^po,oE 2 ) 



+ iW(p 3)4 S 3 + PlfiES) + f;NN(p2 A S 2 + P0,2E) 
+ f;NNN Pl,±S + f;NNNN P0,4 ■ 

A careful application of the algorithm (A.3)-(A.5) yields therefore 



Po,o = 



P2,2 



P4,4 



1 1 

"2 + vT+f T 



n 



Pl,2 



P3,4 



r 2 



(i + r 2 ) 3 / 2 ' 

1 y^F 

2 (1 + r 2 ) 5 ^ 
i 

i - 



2 - 



VT+r 2 ' 
(2 + 3r 2 ; 



P0,2 = 



P2,4 = 

Pl,4 = 



3r 4 f (i + r 2 ) 3 / 2 

| + ^(VTTp-i) 
i i 



rilVi + r 2 -! 



2r 2 v / TTr 2 



(A.6) 



(A.7) 

(A.8) 
(A.9) 
(A.10) 
(A.ll) 
(A.12) 
(A.13) 
(A.14) 
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Po ' 4 64 



(A.15) 



Note that all these new universal functions have a smooth limit at V = 0, where they 
coincide with the Robin coefficients: ^pK, — respectively. 

These new universal functions can be computed by hand after noticing that, by virtue 
of (A.3)-(A.5), one deals repeatedly with the integral 



oc 



R{a^Y) = du ■ i n . r2 w (A.16) 

Jo [u 2 + u + 4(1 + T 2 )\ b 

For example, one finds 

i?(4, 7/2, r) = R(2, 5/2, T) - i?(3, 7/2, r) - 1(1 + T 2 )R(2, 7/2, T), (A.17) 
and this leads to the result (A.15), because 

R(2, 5 /2,r) = -l^ + l^, (A,8) 

i?(2,7/2,r) = _±Atf(2,5/2,r), (A.19) 
R(3, 7/2, T) = R(l, 5/2, T) - R(2, 7/2, r)- 1(1 + r 2 )i?(l, 7/2, T). (A.20) 
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